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We develop a theory which describes the behaviour 
of eigenvalues of a class of one-dimensional random non- 
Hermitian operators introduced recently by Hatano and Nel- 
son. Under general assumptions on random parameters we 
prove that the eigenvalues are distributed along a curve in the 
complex plane. An equation for the curve is derived and the 
density of complex eigenvalues is found in terms of spectral 
characteristics of a "reference" hermitian disordered system. 
Coexistence of the real and complex parts in the spectrum 
and other generic properties of the eigenvalue distribution for 
the non-Hermitian problem are discussed. 

PACS numbers: 72.15Rn, 74.60Ge, 05.45. +b 



1. Complex eigenvalues of non-Hermitian random 
Hamiltonians have recently attracted much interest 
across several areas of physics However, the actual 

progress in understanding the statistics of such eigenval- 
ues is mostly limited to the so-called 'zero-dimensional' 
case, i.e. to random matrices with no spatial structure. 
In this case a number of models can be treated analyt- 
ically and their basic features are relatively well under- 
stood, although on different levels of rigour (see, for in- 
stance, recent works [^-^] and references therein). In 
contrast, little is known about spectra of non-Hermitian 
Hamiltonians in one or more dimensions. One of the chal- 
lenging problems here involves an unusual localisation- 
delocalisation transition predicted by Hatano and Nelson 

0. 

Motivated by the studies of statistical mechanics of the 
magnetic flux lines in superconductors with columnar de- 
fects, Hatano and Nelson considered a model described 
by a random Schrodinger operator with a constant imag- 
inary vector potential. Appealing to a qualitative rea- 
soning they argued that already in one dimension some 
localised states undergo a delocalisation transition when 
the magnitude of potential increases. The eigenvalues 
corresponding to the localised states are real and those 
corresponding to the extended states are complex. The 
results of numerical calculations presented in |jj support 
these conclusions. They also show a surprising feature of 
the eigenvalue distribution in the model: the eigenvalues 
are attracted to a curve in the complex plane. 

In this short communication we explain this feature. 
We also derive an equation for the curve and obtain the 
density of complex eigenvalues in terms of the well known 
objects of the conventional (Hermitian) disordered sys- 
tems. 



Most of our discussion will involve the lattice case 
which is technically simpler. Following |l[, we will con- 
sider a one-dimensional non-Hermitian Anderson model 
whose eigenvalue equation reads as follows 
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The relevance of this non-Hermitian Anderson model to 
physics of vortex lines in superconductors is explained 
in pi. In these works two particular cases of the model 
are discussed: (i) 'site-random' model with random q^ 
and with constant hopping elements given by r/k = h and 
£/. = —h [h is real); and (ii) 'random-hopping' model 
with qk = and with random hopping elements given by 
rjk = g + cik and £/. = g~ak, where g ± are some constants 
and the otk are random variables. 

Our basic assumptions about the coefficients in Eq. 
(|]) are as follows: {(qk,£,k,Vk)} is a stationary sequence 
of random three-dimensional vectors such that (log(l + 
l<Zo|))j (Co), {Vo) are finite. The angle brackets denote 
averaging over the disorder. 

2. We start our analysis with a standard transforma- 
tion which is often used in the theory of differential and 
difference equations. Let us put tpk — Wkipk in Eq. (g) 
and choose the weight u>k so that to make the resulting 
equation symmetric. For instance, if we set 



wq = 1 and Wk = e 



— p 2 J2i=0^ j ' 



■V.i) 



if k > 1, 



(3) 



then this transformation reduces the eigenvalue problem 
(|lf)-(|2|) to the following one (we will use the notation Ck = 
exp[(£fc + 7?fe)/2] onwards) 



Ck-i<Pk-i - Ckifk+i + qk^k = z(p k 



<J n+1 Wl(pi, (fin = W n lf . 



(4) 
(5) 



From now on, we will deal with Eqs. (f|)-(|5|). [Obviously, 
the eigenvalues of (|l])-(||) an d coincide.] 

We can now introduce a "reference" symmetric eigen- 
value problem which will be used in our analysis. This 
reference problem is specified by the same Eq. (g) and 
with the following boundary conditions (b.c.) 



fn+l =0, lfi = 0. 



(G) 



One can rewrite the eigenvalue problems (Q),(|^) and 
(g),(|^) in the matrix form Hip — zip and Hip = zip, 
respectively. We use the calligraphic H. for the non- 
symmetric problem (||),(||) and the usual italic H for the 
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symmetric problem (Q) , (||) . H is a symmetric tridiagonal 
nxn matrix with the {qk} on the main diagonal and the 
{— Cfc} on the sub-diagonals. H is "almost" tridiagonal: 
the only non-zero elements of the difference V = H — H 



are V ln 



u n e io and V nl = 



-wiw n 1 e ri ™ 



Let zi, . . . ,z n denote the eigenvalues of Ti.. Distribu- 
tion of {zj} in the complex plane z = x + iy is described 
by the measure 

1 ™ 

dv n (z) = dv n (x,y) = - J~]5(x- Xj)S(y - yj)dxdy. 
U 3=1 

Our goal now is to find the limit of dv n (z) when n — ► oo. 
To do this we will calculate the "electrostatic" potential 
F(z) = F(x, y) of the limit distribution of the eigenvalues 
of 7i: 



F(z) = lim F n (z), 



(7) 



where 



F n {z) 



log \z - z'\dv n {z') = - log | det(W - zl)\ (8) 



is the potential of dv n (z). In Eq. (ph / is the nxn identity 
matrix. 

According to the potential theory, the existence of 
the limit in Eq. (Q) for almost all z implies the exis- 
tence of the weak limit dv{z) = linin^oo dv n (z). The 
limit measure is determined by the Poisson equation, i.e. 
du{z) = dv{x,y) = l/(2w)AF(x,y) dxdy §|. 

The idea of using potentials to study eigenvalue distri- 
butions is not new and goes back to the 1960s at least, to 
studies of the eigenvalue distribution of Toplitz matrices 
H . In the context of random matrices this idea has been 
used since works ||[l(|. 

We will prove below that with probability one the limit 
in Eq. (0) exists for almost all z and is given by 



F(z) 



a if <f>(z) < a 
$(z) if $(z) > a, 



(9) 



where a = max((£o), (vo)) an d is the potential of 

the limit distribution of the eigenvalues of H . In other 
words, $(z) = linin^oo $„(z), where 

$ n (z) = ~\og\det(H-zI)\= f + \og\z - \\dN n (\) 
n J-oo 

and N n (X) is the number of eigenvalues of H in the inter- 
val (— oo, A) divided by n. It is well known [ni]Jl3] that 
there exists a continuous non-random function N[X) such 
that with probability one linin^oo N n (\) = N(X) and 
hence 



log|A-z|diV(A). 



(10) 



We will call the set which supports dv(z) the limit spec- 
trum of Ti.. $>(z) as a function of x = Rez and y = Imz 



is harmonic everywhere in the complex plane except that 
part of the real axis which supports dN(\). Therefore, 
Eq. (||) says that the complex part of the limit spectrum 
is determined by the equation $(z) = a. This equation 
defines a curve in the complex plane. We denote this 
curve by C. The density of the eigenvalue distribution 
with respect to the arc-length measure ds on this curve 
equals to (2tt)~ 1 times jump in the normal derivative of 
F across the curve. Computing the derivative gives 



dv 
ds 



1 

27 



-foo 



dN{\) 



G C. 



(11) 



N(X) is the classical integrated density of states as- 
sociated with Eq. (Q). It is very important that Q(z) 
coincides, up to an additive constant, with the Lya- 
punov exponent j(z) of the same equation. Namely, 
the well known Thouless formula jl4| states that $>(z) — 
-f(z) + |((£o) + (vo)) (for a rigorous proof of this formula 
and further references see e.g. EH). 

The spectral properties of Ti. are, to a large extent, de- 
termined by the behaviour of 7(2). This will be discussed 
later. Here we only mention that the equation $(z) = a 
which defines the curve C is equivalent to 



7(*) = ^K£o>-<t»)|. 



(12) 



Eq. ( fl2| ) is more transparent. Loosely speaking, it says 
that a non-real z belongs to the limit spectrum of Ti. only 
if there exists a solution tp of Eq. (|j) such that \(p n \ and 
w~ l are of the same order when n — > 00 [compare this 
with Eq. (||)]. As an immediate consequence of Eq. ( |T^ ) 
we obtain the following. If (£0) = (i]o) then there is no 
complex part in the limit spectrum of 7i, i.e. dv(z) is 
supported on the real axis only. 

To derive Eq. (^|) let us consider a non-real z such that 
$(z) 7^ a. Then G = (H - zl)^ 1 is well defined and 
H-zI={I+ VG){H - zi). Therefore 

F n (z) = $ n (z) + -log\d n (z)\, 
n 

where d n (z) = det(I + VG). Expanding this determinant 
yields 

d n (z) = [1 + V ln G nl ][l + V nl G ln ] - V ln V nl G n G nn , 

where Gjk denotes the (j, k) matrix entry of G. To prove 
Eq. @ it remains to find the linin-.oo n _1 log |d„(2)|. 
Using Gi„ — G n i — n?=i c k/ det(Ho-zI) one obtains 

|Vi„G„i| = exp{n[(£ ) - $(2) + r n {z) + s n ]}, 
\V nl G ln \ = exp{n[(?7o) - $(z) + r n (z) + t n ]}, 

where r n (z) = $(z) - $„(z), s„ = rc -1 X)£=o & - (£0), 
and t n = n > z2™=iVj — (^o)- As mentioned above, 
*&n(z) ~~ > $(2) and hence r„(z) — > when n — * 00. 
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Also, because of the ergodic theorem lim r . 
lim^oo t n = 0. 

Consider first the case when 3>(z) < a. Without loss 
of generality we may suppose that a = (?/o)- Then 
VniGin grows exponentially with n. On the other hand, 
$(z) > (£ ) for every non-real z |ll| and Vi n G n i tends to 
zero exponentially fast when n — > oo. Besides, Vi n V n i — 
exp[(^o + ?7n)/2] and the stationarity of the sequence {r] n } 
implies that 



A k = 



lira n- 1 log\V ln V nl \=0. 



(13) 

$(z), hence 



Therefore lim„_ ) . o0 n 1 \og\d n (z)\ = 
linin^oo F n (z) = a if 3>(z) < a. 

Consider now the case when $(2;) > a. Then both 
V n \G\ n and Vi n G n i are exponentially small if n is 
large. Therefore the "dangerous" term in d n (z) is 1 — 
VxnVniGnGnn- It can be shown |16j that a < argGn < 
it — a, where a is strictly positive (0 < a < tt/2) and de- 
pends only on z and on the realisation of {(<ft, ijk)}- 
One can deduce from this inequality that 

|1 - Vi n V n iGiiG nn \ > sina. 

This estimate together with Eq. ([l3]) and the simple in- 
equality \Gjj\ < llmz^ 1 prove that log |d„(z)| — > 
when n — ► 00. Thus lim,^,^ F n (z) — &(z) if &(z) > a. 

Since the limit distribution of eigenvalues of TL is one- 
dimensional, one can also use the conventional tech- 
nique of Green's functions to obtain this limit distribu- 
tion. Let z be non-real and f n (z) = n^ 1 tr(7i — zT) , 
g n (z) = n _1 tr(iJ- f(z) = lim^oo /„ (z) , and 

g{z) — limn^oo g n {z). Revising the above argument one 
can show that f(z) = if <J>(z) < a and f(z) = g(z) if 
$(z) > a. This can also be shown by working directly 
with the Green's functions f n (z) and g n (z). 

If {(<?&, r/k)} is a sequence of independent identically 
distributed random vectors with finite variance and 
rjk are not necessarily independent), then more precise 
information on the large-n behaviour of r n {z), s n , and t n 
is available. Revising the above analysis and employing 
the limit-theorem techniques one can show in this case 
that 



max dist(zj, £) 

Zj :| Im Zj \>e 



(14) 



is of the order 1/y/n for any positive e. 

3. In this section, we consider another approach to 
the eigenvalue problem (0)-@. It is based on the theory 
of Lyapunov exponents of products of random matrices. 
This approach explains the true reason for the appear- 
ance of the complex eigenvalues and gives an independent 
derivation of Eq. (|l2|) . It also provides more precise in- 
formation about the finite-n behaviour of the eigenvalues 
of H. Finally, it easily extends to the case of differential 
equations. It does not provide though the explicit form 
of the limit distribution derived above. 

Let us introduce the usual transfer-matrices 



f-'A: 



Ck 







(15) 



Then the solution of Eq. (|j) with initial data (ipo,(pi) 
can be written as 



(16) 



where Sk(z) = AkAk-i ...A\. Thus the eigenvalue prob- 
lem (0) - (0) reduces to 



(B n S n (z) -w^l) (Lp 1 ,(f ) T = 0, 



(17) 



where / is the 2x2 identity matrix and B n is a 2 x 2 diag- 
onal matrix, B n = diag{exp[-i(£ - rjo)], exp[-|(£„ - 
rj n ))}. In other words, the eigenvalues of TL solve the 
equation det(B n S n (z) — w~ +l I) = 0. This equation is 
equivalent to d n (z) = 0. 

The Lyapunov exponent j(z) of S n {z) is defined as 



7(2) = lim -\og\\S n {z) 

n— >oo n 



(18) 



It is well known ( |I7j| , [|11[ ) that under our basic as- 
sumptions on {(qk, £,k, Vk)} the limit on the right-hand 
side in Eq. (|l8|) exists with probability one. 

For every fixed z the right-hand side of Eq. ( |l8| ) gives 
with probability 1 the same non-random number as the 
Thouless formula and the well known Furstenberg's for- 
mula fl7}| . To proceed, however, we have to study 
7(2;) as a function of z when a (typical) realisation of 
w = {(Qk,S,k,r]k)} is fixed: 7(2) = j(z,uj). Thus, we are 
faced with products of random matrices depending on a 
parameter. The latter were studied in ffgj . 

Let E denote the spectrum the operator defined by the 
left hand side of Eq. (|) on ^(Z) [this set can equivalently 
be described as the support of the measure dN(X)]. It 
was proved in |l8| that with probability 1 the conver- 
gence in ( |l8| ) is uniform in z if \z\ < C and z does not 
belong to an e-neighbourhood of E; here C is an arbi- 
trary fixed large number and e is an arbitrary strictly 
positive number. It can be deduced from this fact that: 
(a) if /J, n {z) is the largest by modulus eigenvalue of B n S n 
then linin^oo n~ l log |/i n (^;)| = 7(2) and the convergence 
here is again uniform in z in the same domain (clearly, 
the other eigenvalue is e^ n+Vo / fi n (z)) . (b) The following 
function 7(2;) = 7(2; -M0) is now defined for almost every 
sequence to simultaneously for all real x, —00 < x < 00; 
the latter limit exists because 7(2) is monotone in y (the 
property which is obvious from the Thouless formula). 

There is a subtlety here. Namely, we can consider 7(2;) 
defined by ( |l8|) for every fixed(!) x but not for all x 
simultaneously while the potential is fixed. Moreover, it 
turns out that if a typical sequence lo is fixed then there 
exists a " large'' subset O w of E such that for all x € O w 
the limit in (|18j ) either does not exist or does not coincide 
with 7(2;) ]18[ |. However the following statement which 
is very important for us holds true: 



7(1) > lim sup -log 1 15„ (i)||=7(a;) (19) 

n— >oo n 
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A more detailed discussion of properties (a), (b), and jl9| ) 
along with a very simple proof of the mentioned above 
uniform convergence in ( |l8| ) will be provided in filjfl . 

Combining property (a) with the statement that 
lim„^oo rT 1 logiUn+i = ((£o) - ^°H/ 2 one conclu des 
that all non-real solutions of Eq. (|17|) are asymptoti- 
cally (as n — > oo) attracted to the curve L given by Eq. 
(|l2|). Moreover, it follows from the uniform convergence 
of u n (z) that linin^oo p n — 0, where p n is defined by 

try 

It is worth noticing that Eq. ( |18| ) can very rarely be 
used for calculation of the Lyapunov exponent "f(z). In 
order to make use of Eq. ( p"2| ) one has to compute 7(2). 
One way to do this is provided by the Thouless formula. 
The other one is the well known Furstenberg's formula 
p7[ . To a large extent, both approaches are equiva- 
lent and in order to use any of these formulae one has 
to find an invariant measure of a certain transformation 
for which purpose an integral equation has to be solved. 
Though in general calculations of this sort are difficult, 
for some distributions the explicit expression for 7(2) and 
N(X) can be found; one of these examples is discussed 
below. 

The properties of the two functions, 7(3;) and 7(3;), are 
crucial for further analysis. Suppose for simplicity that 
in addition to our basic assumptions {q n ,^niVn) form a 
sequence of independent identically distributed random 
vectors. Then it is easy to show that for sufficiently large 



log|x|-d <j(x) <log|x|+Ci, 



(20) 



where C\ depends only on the distribution of (qo, £0, 770)- 
Then, obviously, there exist C such that all solutions of 
Eq. ( |l2| ) belong to a circle \z\ < C. 

The other useful property of 7(3;) is as follows: if in 
addition to the above assumptions qo takes at least two 
different values then 7(12;) > is strictly positive (see 
^7| , ^2|Jl3| ) . Finally, under these conditions, 7(3;) is a 
continuous function of x (see, e.g. ]jl|). 

We are now in a position to describe the curve C. Con- 
sider the inequality 



(21) 



Because of the continuity of 7(3;), the solution of this in- 
equality is given by a union of disjoint intervals: Lij [dj , at-] 
with cij < a'j. Next, for every x £ [aj,a'j], consider a 
positive y = yj(x) such that the pair (x,y) solves Eq. 
(|l2|). This solution exists and is unique because if x 
is fixed then 7(2; + iy) is strictly monotone continuous 
function of positive y and lim y ^ +00 7(3; + iy) = 00. The 
curve C (the solution of Eq. (|l2|)) is a union of discon- 
nected contours: C = DjCj. Each contour Cj consists of 
two symmetric smooth arcs y = yj(x) and y = —yj(x), 
where dj < x < a^-; the points (dj, 0) and (a'j, 0) are 
the end-points of these arcs. [We notice that it is easy 
to construct examples with a prescribed finite number of 



contours. In general, there is no obvious reason for the 
number of contours to be finite for an arbitrary distribu- 
tion of (£0, Vo, Qo)-] 

Our next goal is to describe the behaviour of the real 
eigenvalues of Eqs. (|^)-(|^). It turns out that this be- 
haviour is governed by the following remarkable property 
of 7(3;): this function is upper semi-continuous at each 
point x where 7(3;) = 7(3;) (this equality holds for almost 
every x; remember that 7(3;) depends strongly on u> [fi8|) 

It can be deduced from this property jl6| that for every 
strictly positive e the spectrum of TL lies outside a domain 
surrounded by the following strip: 



Vj. £ = {z e C 



dist(|z|,£j) < e} 



if n is large enough. In other words, with probability 1 
the spectrum of TL is wiped out from the interior of every 



contour Cj as n 



Finally, here is a description of the limit spectrum of 
Eqs. (^)-(||). Let E be the same as above (i.e. the spec- 
trum of the reference symmetric problem) . Then the set 



£u{s\ , («,.«','} 



(22) 



is the limit spectrum of TL. 

4. We start our discussion with a remark on the spec- 
trum of the limit operator TL defined by the left-hand 
side of Eq. ([!]) on ^(Z). For simplicity we suppose that 
all coefficients in Eq. (fj]) are bounded. It turns out that 
for a wide class of distributions of {(£fc, f]k,qk)} the spec- 
trum of TL is a two-dimensional subset in the complex 
plane and the limit spectrum of TL given by Eq. ( p2| ) is 
embedded into this set. This phenomenon seems to be 
surprising because TLip converges to TL(p when n — > 00 (ri 
is the dimension of TL) for every cp S ^(Z). 

Since the problem is non-Hermitian, the spectrum may 
depend on the choice of boundary conditions (b.c). In- 
deed, all the eigenvalues of Eq. (|l|) with the Dirichlet 
b.c, i.e. when ip n +i = ipo = 0, are real. It is remark- 
able, however, that the boundary conditions of the form 
(VVi+i, VVi) T — vpo) T , where B is a fixed real non- 

degenerate 2x2 matrix, lead to the same Eq. ( |l2|) (re- 
gardless of the choice of B) in the limit n — > 00. For 
diagonal B this fact can be readily seen from our proof 
of Eq. (|^) . The general case of a non-degenerate B is less 
transparent and our proof of Eq. (|l2] ) in this case fill 
relies on the above mentioned properties of products of 
random matrices . 

The reader might have noticed that our derivation of 
Eq. (^) is based only on the fact of existence of N(X). 
(Under our basic assumptions the existence of N(X) is 
ensured by the stationarity of the coefficients.) Thus one 
can easily extend our argument to other classes of coeffi- 
cients. 

For instance, instead of a random sequence 
{(£fc> Vk, qk)} one can take a periodic one (i.e. the oj., 
and r\k all have a common period). Then Eq. (M) 
holds with the following obvious change. In the periodic 
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case the angle brackets denote averaging over the period. 
The complex part of the limit spectrum is described by 
the same equation ( |l2| ) in which 7(2) is the Lyapunov 
exponent of the symmetric reference equation (^) whose 
coefficients are now periodic. 

The geometry of the limit spectrum of TL in the peri- 
odic case follows two patterns. If (£) = (77) then the limit 
spectrum is real and coincides with E, the spectrum of 
the symmetric reference equation. If (£) ^ (77) the limit 
spectrum is purely complex, i.e. the limit distribution of 
eigenvalues is supported on the curve defined by Eq. (12). 
It is worth mentioning that in the periodic case this curve 
is a union of a finite number of analytic contours. Indeed, 
because 7(2) = on £ in the periodic case, the above 
mentioned arcs join up smoothly. In either case (real or 
complex spectrum) the corresponding eigenfunctions are 
extended. The case when rjk, qk)} is quasi-periodic 
is much more delicate and deserves a separate study. 

The rest of our discussion is based on known results 
on Hermitian disordered systems and the formulae de- 
rived above. From now on we assume that the {qk} is a 
sequence of independent identically distributed random 
variables which is also independent of the {c^}. It is 
instructive to consider first an exactly solvable model. 

Let Cfc = 1 and qk are Cauchy distributed, i.e. 
Prob{qfc G A} = 7T _1 J A dq b/(q 2 + b 2 ). In this special 
case an explicit expression for -/V(A) is known (see, e.g. 
IP^I ). Using this expression and the Thoulcss formula 
one obtains that 



4cosh 7 (z) = y/(x + 2) 2 + (b+\y\) 2 + 
^(x-2) 2 + (b+\y\) 2 , 



(23) 



where x = Re z and y = Im z. One can use this result to 
obtain the limit spectrum of a family of non-Hcrmitian 
operators TL (fi_])-(g). Suppose that in Eq. (pj) r]k = — £fc 
for every k (e.g. one can take r\k = h and = —h as in 
the site-random model). Then in the reference equation 
Ck = 1 and by straightforward computations involving 
Eqs. (12) and ( p3|) one finds that the complex part of the 
limit spectrum of TL consists of the two arcs 



y(x) = ± 



(K 2 - 4)(K 2 - x 2 ) 



K 2 



-b 



-Xb < X < Xb, 



(24) 



where K = 2 cosh (770) and Xb is determined by the condi- 
tion y{xb) = 0. The real eigenvalues of TL are distributed 
(in the limit n — > 00) in the intervals (— 00,— Xb) and 
(xb, +00) with density equal to the density of eigenvalues 
of Eq. (0) in these intervals. 

The 6-dependence of the complex part of the limit spec- 
trum is remarkably simple. If b — the arcs (|24|) form the 
ellipse x 2 /K 2 + y 2 j{K 2 — 4) = 1. As b increases each of 
the arcs moves (by translation in the y-direction) towards 
the real axis and reduces in length. At b = \/K 2 — 4 the 
arcs (hence the complex part of the limit spectrum) dis- 
appear. In other words, if K < K cr = V^ + ¥ Eq. (EI 



cannot be solved and the limit spectrum of TL is real. If 
K > K cr then the limit spectrum of TL has always two 
branching points from which the complex branches grow. 

We should mention another example of the reference 
eigenvalue problem (^) amenable to analytic treatment. 
In this example the Ck are independent exponentially dis- 
tributed random variables and qk = 0. Thus this example 
is related to the random-hopping non-Hermitian Ander- 
son model. The known expression for 7(2) |l9| ] is not 
that simple as in the previous case and in order to anal- 
yse Eq. ( p"2| ) one has to employ certain approximations. 
This analysis goes beyond the scope of the present pub- 
lication. 

The next part of our discussion involves general anal- 
ysis of Eqs. (H) and (|l2|). Denote 2g = (770) — (Co)- As 
mentioned above the limit spectrum of TL is entirely real 
if g = 0. First, we want to show that in the case when 
qk = (random hopping model) the limit spectrum has 
a non-trivial complex part for every non-zero g. Indeed, 
in this case the transfer matrices Ak (jlj) have zero di- 
agonal elements when z = 0. Hence their products and 
the norm of S n (0) can be easily computed (see, e.g. Q). 
From this computation and the stationarity of the c& it 
follows easily that 7(0) = 0. Therefore, if g ^ Q the solu- 
tions to the equation 7(2;) = \g\ constitute a non-trivial 
(non-real) curve. 

This is not the case if the distribution of go (hence of all 
qk) is non-degenerate (i.e. go takes at least two different 
values). Then there exists g C r > such that for all 
1 5 1 < g^r the limit spectrum of TL has no complex part. 
Indeed, if the qk are non-degenerate 7(2) = j(x + iO) is 
strictly positive on the limit spectrum E of the reference 
eigenvalue problem (Q),(S). The continuity of 7(2), Eq. 
( P0| ) and the inequality ~{x + iy) > j(x) imply that 

9& = min7(x) > 0. 

If in addition the coefficients in Eq. (Bj) are bounded, i.e. 
cl + qf < C for all k, then E is a bounded set. Therefore 

9 { cr = 

is finite. If \g\ > gi? the inequality ( pl| ) is satisfied for 
every point of E. Hence the limit spectrum is purely com- 
plex in this case. If either cq or qo takes arbitrary large 
values with non-zero probability, then E is an unbounded 

,( 2 ) 



set and in view of Eq. fl20|) g cr = +00. 

In summary, if \g\ < g^ the limit spectrum of TL is 

(2) 

entirely real, if \g\ > g cr the limit spectrum is entirely 

complex and if gcV < \g\ < gh) the limit spectrum has 
real and complex parts. In the latter case the branching 
points from which the complex branches grow out of the 
real eigenvalues are determined by 7(2;) = \g\. 

The limit distribution of the real eigenvalues of TL is 
described by N(X), the integrated density of states of 
the reference symmetric equation. The limit distribu- 
tion of the non-real eigenvalues of TL is described by Eq. 
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( |ll| ) . It should be noticed that the density of the non-real [20] R.A.Janik, 



eige lvalues given by Eq. ( 11 ) is analytic everywhere on 



mat/9705098 



M.A.Nowak 
P.W 



C except the (real) end-points of the arcs. (If the limit 
spectrum is entirely complex then this density is analytic 
everywhere.) The behaviour of the density of the non- 
real eigenvalues near an end-point Xj of an arc depends 
on the regularity properties of N(X) at this point. If 
the density of states AT' (A) of the reference equation is 
smooth in a neighbourhood of the point A = Xj then the 
density of the limit distribution of the non-real eigenval- 
ues of Ti has a finite limit as z approaches Xj along the 
arc. Also, in this case the tangent to the arc exists at Xj 
and is not vertical. In other words, if N'(X) is smooth in a 
neighbourhood of a branching point A = Xj the complex 
branches of the limit spectrum grow out of Xj linearly. 
This may not be the case if AT' (A) is not smooth. 
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